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ABSTRACT 

Using the observed submillimetre source counts, from 250-1200 microns (including the 
most recent 250, 350 and 500 micron counts from BLAST), we present a model capable 
of reproducing these results, which is used as a basis to make predictions for upcom- 
ing surveys with the SPIRE instrument aboard the Herschel Space Observatory. The 
model successfully fits both the integral and differential source counts of submillime- 
tre galaxies in all wavebands, predicting that while ultra-luminous infrared galaxies 
dominate at the brightest flux densities, the bulk of the infrared background is due 
to the less luminous infrared galaxy population. The model also predicts confusion 
limits and contributions to the cosmic infrared background that are consistent with 
the BLAST results. Applying this to SPIRE gives predicted source confusion limits 
of 19.4, 20.5 and 16.1mJy in the 250, 350 and 500 micron bands respectively. This 
means the SPIRE surveys should achieve sensitivities 1.5 times deeper than BLAST, 
revealing a fainter population of infrared-luminous galaxies, and detecting approxi- 
mately 2600, 1300, and 700 sources per square degree in the SPIRE bands (with one 
in three sources expected to be a high rcdshift ultra-luminous source at 500 microns) . 
The model number redshift distributions predict a bimodal distribution of local quies- 
cent galaxies and a high redshift peak corresponding to strongly evolving star-forming 
galaxies. It suggests the very deepest surveys with Herschcl-SPIRE ought to sample 
the source population responsible for the bulk of the infrared background. 
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1 INTRODUCTION 

Over two decades ago, a large number of galaxies that 
emit the bulk of their luminosity in the restframe far-IR 
were detected in the IRAS All-Sky Survey (typically at 
z < 0.1). These so-called luminous and ultra-luminous in- 
frared galaxies (LIRGs IO^Lq < LiR(8-1000/im)< IO^^Lq 
and ULIRGs Lir > 10^^ Lq) are powered by a com- 
bination of star formation and active galactic nucleus 
jSoifer et all Il98i1 ). but only recently have they been 
shown to be an important population in the early Uni- 
verse. This is in part due to the achievements of sub- 
millimetre continuum observations using ground-based fa- 
cilities: pioneering surveys at 850/im with SCUBA on 
the JCMT begat the discovery of submillimetre galax- 



ies (SMGs: e.g.. ISmail. Ivison fc Blain I[l997l. iBarger et al. I 

39|) 



Il998l . iHughes et al. 1 1 19981 . lEales et al. 



19991) which were 



subsequently constrained to be main ly distant star-forming 
galaxies (e.g.. [Chapman et al. |[2003l ). These characteristics 
were shared with SMG s found in other submillimetre bands 
e.g., 1100 & 1200 »m jLaurent et al ll2005l . iBertoldi et all 
|2000D and 350 fim (iKhan et al. Il2007t). Larger su rveys (e.g. 
the SCUBA SHADES survey. iMortier et al. 1120051) have con- 
firme d these sources are strongly evolving (jCoppin et al. I 
I2OO6I ). However, the discovery of SMGs still poses chal- 
lenges to se mi-analytical hierarchica l models of galaxy for - 
mation (e.g. iGuiderdoni et al.lll998l . iBalland et al~ll2003l ). 
and questions remain over their role in the f ormation of ellip- 
tical galaxies and supermassive black holes iMagorrian et al.l 
((1998 ) and the energy budget between star-formation and 
accretion in the Universe. 
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In this work we present a galaxy evolution model that 
successfully reproduces the source counts from 250-1200^m, 
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including both the large area SCUBA surveys and the lat - 
est results from the BLAST telescope (jPascale et al. IboOSh . 
In Section [2] we describe the model and present fits to the 
galaxy counts in Section [S] The launch of SPIRE on-board 
the Herschel Space Observatory offers an opportunity to ex- 
amine an SMG population that overlaps with ground-based 
observations and IR-luminous galaxies selected at mid - far- 
IR wavelengths (e.g., with IRAS, AKARI, Spitzer). SPIRE 
will perform surveys at 250, 350, 500/im and in Section |3] 
we discuss the prospects for upcoming surveys with Her- 
schel. Throughout this work a concordance cosmology of 
Ho = 72kms"^Mpc"\ r2 = 0.3, A = 0.7 is assumed. 



2 THE GALAXY EVOLUTION MODEL 

To model the submillimetre source counts we use a far- 
IR backward e v olution framework fol lowing the models of 
iPearsonI l|2005l '). iPearson et"al~l (120071 ). These models were 
previously successfully used to reproduce the combined 
mid-infrared source counts from ISO & Spitzer at 15/im 
& 24/im. These models have now been updated to pro- 
duce source counts f r om 1- 1000 /im and will be reported 
in detail in iPearsonI (|2009l 'l. Althou gh submillimetre lu- 
mino sity functions are available (e.g. ISerieant fc HarrisonI 
I2OO5I ). to model the counts we retain the 60/im luminos- 
ity fu nction derived from t he IRAS Point Source Cata- 
logue l|Saunders et al. I|2000D since it is defined around the 
peak of the dust emission and contains a large ensemble of 
sources segregated by population class. The source counts 
are fit to the wavelength where the luminosity function is 
defined, \lf, which sets the baseline normalization of all 
parameters. To predict the counts at other wavelengths, 
the luminosity function is shifted to the observation wave- 
length, Xobs, using the ratio L[\obs) / L[\lf) , obtained via 
model template spectra, no other priori is assumed. Spec- 
tral templates are drawn from four source populations, com- 
prising normal quiescent galaxies and three star-forming 
groups consisting of, with increasing luminosity, starburst 
galaxies [Lir < 10^^ Lq) , LIRGs, and ULIRGs (mod- 
elled on the archetype Arp220). An additional ACN com- 
ponent (based on the emission from a d ust torus) is als o 
included within the model framework of IPearsonI ([200^ , 
however it is found that AGN do not contribute signifi- 
cantly to the source counts in the submillimetre and al- 
though included, their contribution is not considered in this 
work. The normal galaxy s pectral templates are from the li - 
braries described in Efstat hiou fc Rowan- Robinson I (|2003l ) 
which exhibit cold far-IR/submillimetre colours, with spec- 
tra peaking between 100-200^m. The adopted starburst, 
LIRG & ULIRG spectral template s are taken from the 
spectral models of lEfstathiou et al. I (|2000 h. which provide 
go od fits to the IRAS, ISO and S p itzer galaxy p opulations 
( [ Rowan-Robinson et alTI |2004| . iRowan- Robinson et al. I 
I2OO5I ). Note that all templates are independent of the ob- 
served data sets being fitted. 

Follow-up SCUBA imaging of local IRAS-selected 
galaxies has implied colder far-IR-submillimetre colours in 
SMGs than those derived from galaxy spectra based purel y 
on IRAS colours (|Dunne et al. ll2000l . [Vlahakis et al. Il2005h . 
The colours of our model templates agree with this, as they 
follow the trend of the local galaxy colours extremely well 
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Figure 1. Colour-colour distributions of the model normal, star- 
burst and U/LIRG templates compared with the local submil- 
limetre & IRAS far-infrared colours from Dunne et al. (2000), 
Vlahakis et al. (2005), with locus of normal (DISK), starburst 
(SB) and the ULIRG ARP 220. The markers along the SED tracks 
correspond to redshift steps of (5z=0.1. 



in Figure [T] Although deeper SCUB A surveys ar e expected 
to principally select LIRG/ULIRGs ijBlain et al. t,200 2). this 
local sample also comprises lower luminosity starburst and 
cooler no rmal galaxies (also predicted to contrib ute at higher 
redshifts lEfstathiou fc Rowan-Robinson II2003I ) . 

The star-f orming populat i ons follow the burst evolution- 
ary scenario of lPearso 3 (|2005^ ■ |Pearson (.20Q9 ), modelled by 
an exponential function to « ~ 1 and a power-law there- 
after. This evolution is consistent with a rapid onset of star- 
formation at high redshift, a gradual decline to redshift of 
~1 and a sharp decline in activity to the present epoch. 
The relative contribution of each component to the overall 
star-formation rate follows a downsizing pattern with red- 
shift in which the most massive galaxies formed stars at an 
early epoch, thus domi nating the star-formati on history in 
the early Universe fe.g. iMobasher et al. |[2009l ). 



3 SOURCE COUNTS AND RESULTS 

The model fits to the observed source counts at 250, 350, 
500, 850 and llOO^m are shown in Figure [2] with the total 
model source counts shown alongside the respective contri- 
butions of the assumed galactic populations (normal and 
starburst galaxies, LIRGs and ULIRGs). 

Figure [2] panels a, c Ez d show the model fits to the dif- 
ferential counts (normalised to a Euclidean universe) from 
the Ballo on-borne Large- Aper ture Submillimeter Telescope 
(BLAST. iPascale et al. II2OO8I ) survey in the GOODS field 
dPevlin et al. II2009I ). for the 250, 350 & 500/im bands re- 
spectively. In all the BLAST bands, it is predicted that 
the brightest counts (>lJy in the non-evolving Euclidean 
regime of the counts) will be dominated by quiescent normal 
galaxies, expected to be bright, local galaxies at redshifts 
<0.5 (see Figure |3}. Here the model predicts 2.7 sources 
at 250 /xm over the BLAST survey area of 8.7 deg^, com- 



pared with the th ree sources found in the brightest bin of 
the source list of iDevlin et al. I (|2009l ). The steep depar- 
ture from Euclidean counts is caused by the ULIRGs but 
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Figure 2. Model fits to the observed Submillimetre counts from 250 - 1200^m. The total contribution and individual components 
corresponding to the normal, starburst, LIRG and ULIRG populations are overplotted. (a,c,d) Model fits to the observed BLAST 
counts at 250, 350, 500/im from Devlin et al. (2009). Also shown is the SHARC-II 350/im observation from Khan et al. (2007). Source 
counts are differential normalized to a Euclidean universe. (b)The Observed SSO^mi integral source counts from the SHARC-II survey 
of the Bootes field by Khan et al. (2007); also plotted are the normalized 450|am SCUBA source counts from Small et al. (2002). (e) 
Observed 850^m integral source counts from the various surveys by SCUBA (Small et al. (1997), Hughes et al. (1998), Bales et al. 
(1999), Barger et al. (1998), Blain et al. (1999), Small et al. (2002), Cowie et al. (2002), Scott et al. (2002), Knudsen et al. (2006)); and 
the largest SCUBA survey, SHADES, covering ~0.25 sq.deg. Coppin et al. (2006). (f) The SHADES differential counts normalized to a 
Euclidean universe, (g) Observed 1100 [am integral source counts are from the BOLOCAM instrument Black shaded area from Laurent 
et al. (2005) and the MAMBO survey of Greve et al. (2004) normalized from 1200 |am. (g) Differential source counts normalized to a 
Euclidean universe for the AzTEC observations of Perera et al. (2008). 



at the peak of the differential source counts the less lumi- 
nous LIRGs are the dominant population. At 250 /im, the 
model slightly over-predicts the source counts at >60 mjy, 
but this is within the BLAST error bars. There is a sharp 
rise in the counts at the 200 mJy level, and a turn-over be- 
tween 100 and 20 mJy (although the BLAST counts in this 
region may be less reliable as the instrument is confusion- 
limited), with the model predicting a second turn-over at 
fainter flux densities (<10mJy). Due t o the strong nega- 
tive K -corrections in the submillimetre (jFranceschini et al. I 
[l993), the flux densities of distant galaxies are enhanced 
such that the luminosity function at lower luminosities is 
sampled at fainter flux densities, with any break in the 
counts being attributed to a change in the dominant popu- 
lation. The BLAST counts are derived from a P(D) analysis 
rather than source catalogues and provide a statistical con- 
straint on the slope of the source counts at faint fluxes which 
are already source confused. Encouragingly at 350/.im, the 
faintest BLAST counts are consistent with the differential 
counts from the deeper (non-conf used) 350 Atm s urvey using 
SHARC II in the Bootes field by iKhan et all (j2007l ). Fig- 
ure [3) shows the 350 /xm integral source counts from same 
survey and the SCUBA 450Aim counts (|Smail et al. II2002I : 
assuming an Arp 220 spectral template to transform the 
counts to this band). The model fits these observations well, 
predicting breaks in the source counts at ~40 and ~10 mJy, 



and that the deeper SHARC II results are dominated by 
LIRGs. In the 500/im BLAST band, the model fit is excep- 
tionally good, from the steep rise from Euclidean values at 
S <300 mJy, to the turn-over between 30-6 mJy. The model 
predicts turnovers in the counts at fainter fiux densities of 
10, 8, 5mJy in the 250, 350 & 500/im bands, all within 
the constraints imposed by equating the integrated surface 
brightness of the BLAST sources to the emission from the 
infrared background derived from a power-l aw extrapola- 
tion a nd naive cut-off of sources estimated by iDevlin et al. I 
(j2009h to be 7.0±1.3, 7.2±1.7 & 4.6±1.2mJy at 250, 350 & 
500/im respectively. 

At longer submillimetre wavelengths, the model fits are 
compared with the observed integral source counts from the 
myriad surveys carried out with SCUBA at 850 (Figure 
[21;). These observations span two orders of magnitude in flux 
density and thus provide the best pre-Herschel constraints 
on the galaxy counts. The models provide a good flt to the 
counts from the brightest flux densities down to 0.5 mJy - 
below the SCUB A-850 nm confusion limit of 2mJy (from the 
lense d surveys of lSmail. Ivison fc Blain Ill997l . ISmail et al. I 
I2OO2I ). At these levels, due to the strong negative K- 
corrections, we expect to be able to observe relatively moder- 
ate starburst galaxies. At brighter flux densities, ~10mJy, 
ULIRGs are the dominant population (although a signifi- 
cant increasing contribution from normal galaxies cannot be 
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ruled out (jEfstathiou fc Rowan- Robinson 1 12003| ). but from 
the model they are predicted to dominate at >~50mJy). 
The la rgest 850 ^^r^^ survey t o-date (the ~0.25deg^ SHADES 
survey iMortier et al. I boosi l detected 120 sources, effectively 
doubhng the number o f known SMGs The SHADES dif- 
ferential source counts (|Coppin et al. I [20061 ) are shown in 
Figure [2] (f). The best-fitting model requires a break at ~4- 
6 mjy and it is difficult to simultaneously reconcile this with 
the bright-end c ounts, using even the most rec ent evolution- 
ary models (e.g. iRowan- Robinson et al. II2OO9I ). However our 
model fits both the bright-end counts and this break due to 
the inclusion of the intermediate LIRG population, which 
are often omitted in contemporary source count models, be- 
tween the starburst and ULIRG populations. The break is 
predicted to be due to the emergence of these strongly evolv- 
ing galaxies, with their contribution peaking at ~2mJy in 
the differential counts. 

In Figure the integral source counts at millimetre 
wavelengths for the surveys with the BOLOCAM instru- 
m ent at llOO^tm (from the maximum likelihood analysis 
of iLaurent et al. I l|2005l )') and the MAMBO instrument at 
1200/im (normalising the counts to 1100/xm). BOLOCAM, 
MAMBO & SCUBA have aU surveyed the same area - the 
Lockman Hole - and in essence, the counts suggest the mil- 
limetre observations are sampling the same brighter portion 
(Ssso >~8 mJy) of the SCUBA 850 population, expected 
to be dominated by ULIRGs or even HLIRGs (Hyper Lu- 
mmous Infra-Red Galaxies, Lir > 10^^ Lq). This IS sim- 
ply due to the longer wavelengths sampling further down 
the Rayleigh- Jeans slope and therefore preferentially select- 
ing the higher luminosity, high redshift objects. Finally, the 
model fits to t he differential c ounts from the recent AzTEC 
observations of lPereral (|2008l ) in the GOODS fields are pre- 
sented in Figure [2fi. There is a good fit to these observed 
counts, with the higher luminosity sources providing the 
main population. The flattening seen in both the integral 
and differential source counts at fluxes of ~3 mJy is also re- 
produced by the model, representing a shift in the dominant 
population from ULIRGs to LIRGs, with the expectation of 
further flattening below ~lmjy. 



4 THE PROSPECTS FOR HERSCHEL 

The Herschel Space Observatory (jPilbratt II2008I) . launched 
on 14th May 2009, is ESA's next generation infrared mission. 
The Spectral and Photometric Imaging Receiver (SPIRE) 
instrument is one of the focal plane instruments and is de- 
signed for photometr y and spectroscopy between 200-550^m 
ijGriffin et al. I bOOSl ). The three SPIRE bolometer arrays 
(PSW, PMW and PLW, respectively centered on 250, 350 & 
500/im, A/AA ~3, with 139, 88, and 43 pixels) allow simulta- 
neous observations over a FOV of 4' x 8' in the three bands. 
In spire's large map scanning mode the 5a, 1 hour point 
source sensitivities are expected to be 3.7, 5.3 & 4.6 mJy 
for the respective arrays. Despite being near identical to the 
three arrays on BLAST, the 3.5m Herschel primary mirror 
(2m on BLAST) offers superior resolution. 

The ultimate sensitivity of any survey will be the con- 
fusion limit, defined as the threshold of fluctuations in the 
background sky brightness caused by (unresolved) point 
sources below which sources cannot be discretely detected 
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Figure 3. The number redshift distribution for SPIRE bands at 
the confusion Umit of 19.4, 20.5 & 16.1mjy for the PSW 250, 
PMW 350 & PLW 500Atm arrays. Redshift bin size is <5z = 0.1. 



in the telescope beam A/D, where D is the telescope diam- 
eter. The confusion due to faint galaxies is more severe at 
longer wavelengths and smaller apertures and is often char- 
acterized by the number of beams per source, with classical 
limits of 20-40 bea ms per source often adopted r lHoggll200ll . 
IJeong et all2006[ l. The confusion hmits for HerscheLSPIRE 
and BLAST can therefore be compared using the source 
count model: for BLAST the 20 beams per source confusion 
limit is predicted to be 33.7, 33.6 & 23.9mJy in the 250, 350 
fc 500fim ban d s - ag reeing very well with the estimates from 
iDevlin et al.1 (|2009| ) of 33±4, 30±7 & 27±4 mJy (imply- 
ing that the faintest counts reported by BLAST are already 
source confused). For SPIRE, given the larger aperture, our 
models predict 20 beams per source confusion limits of 19.4, 
20.5 & 16.1mJy in the 250, 350 & 500/im bands respectively, 
implying that equivalent surveys with Herschel offer a sen- 
sitivity improvement of 1.5 over BLAST. This advantage is 
significant given the steep nature of the source counts in 
this fiux regime, with the model predicting a shift from the 
ULIRG dominated counts to LIRG domination. Using P(D) 
analysis , the BLAST source counts are able to probe be- 
neath the conventional confusion limit, however P(D) anal- 
ysis is only able to constrain the slope of the source counts. 
The intrinsically deeper images afforded by Herschel should 
be able to resolve and reliably sample this emerging popu- 
lation, which our model predicts will contribute the bulk of 
the Cosmic Infra-Red Background (CIRB). 

Using a similar P(D) analysis, SPIRE surveys could 
also be pushed well below the corresponding SPIRE con- 
fusion limits allowing an actual detection of the turn over 
in the differential counts (note that the BLAST survey can 
only place upper limits). However, surveys with SPIRE to 
these depths are not expected to resolve the sources re- 
sponsible for the break in the source counts (eluded to by 
iKhan et al. 1120071 at 350/im) although this population will 
be accessible using ground-based facilities s u ch as SHARC- 
H/CSO, SCU B A-2/J CMT ,HoUand et al. I (|2006l ). CCAT 
ISebring et "aTl (|2008l '). In the deepest SPIRE confusion- 
limited surveys, we expect 2600, 1300 and 700 sources per 
square degree for the 250, 350 & 500/im bands respectively, 
of which ~12%, 25% & 35% wiU be high redshift ULIRGs. 



Submillimetre surveys: The prospects for Herschel 5 



100 




1 10 100 

S/mJy 



Figure 4. CIRB fraction as a function of flux for the 
SPIRE/BLAST bands, solid-lines are tiie integral and dashed- 
lines are the differential contribution. Also shown are the confu- 
sion limited sensitivities for SPIRE circles and BLAST squares. 



In Figure |3] we show the number redshift distribution for 
the three SPIRE bands at a survey sensitivity correspond- 
ing to the SPIRE confusion limit. In all bands, but most 
predominantly in the short wavelength 250/xm band, a bi- 
modal distribution is seen which can be interpreted as a lo- 
cal contribution from quiescent normal galaxies and a high 
redshift contribution from evolving starburst galaxies, with 
the high redshift peak becoming more pronounced to longer 
wavelengths. The median redshift of the N-z distribution lies 
between 2 < 2: < 3, consiste nt with the redshift dis tribution 
of SCUBA-850Mm sources |Chapman et al. Il2003h . 

Using the model to integrate to the faintest flux levels, 
the total contribution of faint sources to the cosmic infrared 
background in the SPIRE bands is estimated, giving inten- 
sities of 11.0, 6.0 & 2.4nW/mVsr in the 250, 350 & 500^im 
bands. These values agree well with t he COBE -FIRAS re- 
sults of 10.4±2.3, 5.4±1.6 & 2.4±0.6 of lFixen et al. (1991). 
We also estimate that ~80, 85 & 90% of the total back- 
ground resides at z>l at 250, 350 & 500^m respectively. In 
Figure |4] the integral and differential percentage contribu- 
tions as a function of flux for the SPIRE and BLAST bands 
are shown, alongside the corresponding 20-beam confusion 
limits for both instruments. We expect Herschel to resolve 
-30% (70%), 20% (60%) & 12% (45%) at the confusion 
(and optimal instrumental 1-hour integration) level in its 
PSW-250/im, PMW-350 & PLW-500 arrays, with a 
corresponding peak in the background emission to occuring 
at flux densities of 10-25 mjy, 8-20 mjy & 5-10 mjy in the 
respective bands. This implies the very deepest surveys with 
Herschel-SPIRE should sample the dominant source of the 
background, which from our model is expected to be mostly 
LIRGs (rather than more luminous ULIRGs, whose contri- 
bution peaks at slightly brighter fluxes in all wavebands) . At 
faint flux densities (<lmjy) 20-30% of the CIRB remains 
unresolved and the counts in Figure [2] indicate that the 
fainter starburst galaxies will become the dominant popula- 
tion, responsible for the remainder of the total background. 

Therefore we expect the upcoming Herschel SPIRE sur- 
veys to produce the flrst large statistically reliable samples 



of SMGs, taking submillimetre astronomy from the pioneer- 
ing era into one which detailed constraints can be placed on 
the evolution of star-formation in the early Universe. 



5 ACKNOWLEDGEMENTS 

We thank Steve Willner for helpful comments and Andreas 
Efstathiou for providing his galaxy templates. We thank the 
referee for constructive comments that improved this work. 



REFERENCES 

Balland C, Devriendt J. E.G., Silk J., 2003, MNRAS, 343, 107 
Barger A.J.et al., 1998, Nature, 394, 248 
Bertoldi P., et al., 2000, A&A, 360, 92 

Blain A., Kneib J. P., Ivison R.J., Smail I., 1999, ApJ, 512, 87 
Blain A., Smail 1., Ivison R.J., Kneib J.-P., Prayer D.T., 2002, 

Physics Reports, 369, 111 
Chapman S.C., Blain A.W., Ivison R.J., Smail I.R., 2003, Na- 
ture, 422, 695 
Coppin K. et al., 2006, MNRAS, 372, 162 
Cowie L.L.; Barger A.J., Kneib J.-P., 2002, ApJ, 603, 69 
Devlin M. et al., 2009, Nature, 458, 737 

Dunne L., Eales S., Edmunds M.G., Ivison R., Alexander P., 

Clements D. L.,2000, MNRAS, 315, 115 
Eales S. et al., 1999, ApJ, 515, 518 

Efstathiou A., Rowan-Robinson M., Siebenmorgen R., 2000, 
MNRAS, 313, 734 

Efstathiou A., Rowan-Robinson M., 2003, MNRAS, 343, 322 

Fixsen, D.J., Dwek, E., Mather, J.C., Bennett, C.L., Shafer, 
RA., 1998, ApJ,. 508, 123 

Franceschini A., Toffolatti L., Mazzei P., Danese L., De Zotti 
G., 1991, AAS 89, 285 

Guiderdoni, B., Hivon, E., Bouchet, F.R., Maffei B., 1998, MN- 
RAS, 295, 877 

Greve T.R., Ivison R.J., Bertoldi P., Stevens J. A., Dunlop J.S., 

Lutz D., CarilU, C.L., 2004, MNRAS, 354, 779 
Griflin M. et al., 2008, Proc. SPIE 7010, 701006 
Hogg D., 2001, AJ, 121, 1207 
Holland W. et al., 2006, Proc. SPIE 6275, 62751 
Hughes D. et al., 1998, Nature, 457, 616 

Jeong W-S., Pearson CP., Lee, H.M., Pak, S-J.,Nakagawa T., 

2006, MNRAS, 369, 281 
Khan S. et al., 2005, ApJ, 631, 9 
Khan S. et al., 2007, ApJ, 665, 973 
Knudsen K.K. et al., 2006, MNRAS, 368, 487 
Laurent G.T. et al., 2005, ApJ, 623, 742 
Magorrian J. et al., 1998, AJ 115, 2285 
Mobasher B. et al., 2009, ApJ, 690, 1074 
Mortier A.M.J, et al., 2005, MNRAS, 363, 563 
Pascale E. et al., 2008, ApJ, 681, 400 
Pearson CP. 2005, MNRAS, 358, 1417 
Pearson, CP. et al. 2007, Adv.Space Res. 40, 605 
Pearson CP. 2009, in preparation 
Perera E.L. 2008, MNRAS 391, 1227 
Pilbratt G.L., 2008, Proc. SPIE 7010, 701002 
Rowan-Robinson M. et al, 2004, MNRAS, 351, 1290 
Rowan-Robinson M. et al, 2005, ApJ, 129, 1183 
Rowan-Robinson M. et al, 2009, ApJ, 394, 117 
Saunders W. et al. 2000, MNRAS, 317, 55 
Scott S.E. et al., 2002, MNRAS, 331, 838 

Sebring T.A., Radford S., Giovanelli R., Glenn J., Woody D., 

2008, Proc. SPIE, 7012, 70121 
Serjeant S., Harrison D., 2005, MNRAS, 356, 192 
Smail I., Ivison R.J. , Blain A.W., 1997, ApJ, 490, L5 



6 C. Pearson, S. Khan 



Smail I., Ivison R.J., Blain A.W., Kncib J.-P., 2002, MNRAS, 
331, 495 

Soifer, B.T., Neugebauer G., Houck J.R., 1987, ARAA 25, 187 
Vlahakis C, Dunne L., Eales S., 2005, MNRAS, 364, 1253 



